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We present an extensive study of the unique thermal and mechanical dynamics for narrow-line 
cooling on the So - ^Pi ^^Sr transition. For negative detuning, trap dynamics reveal a transition 
from the semiclassical regime to the photon-recoil-dominated quantum regime, yielding an absolute 
minima in the equilibrium temperature below the single-photon recoil limit. For positive detuning, 
the cloud divides into discrete momentum packets whose alignment mimics lattice points on a face- 
centered-cubic crystal. This novel behavior arises from velocity selection and "positive feedback" 
acceleration due to a finite number of photon recoils. Cooling is achieved with blue-detuned light 
around a velocity where gravity balances the radiative force. 

PACS numbers: 32.80.Pj, 32.80.Lg,42.50.Vk, 39.25.+k 



Magneto-optical traps (MOTs) utilizing spin- 
forbidden transitions have recently attracted con- 
siderable attention as starting points for all-optical 
quantum degenerate gases |H, single-system tests of 
Doppler and sub-Doppler cooling theory 2], and es- 
sential components in the next generation of optical 
frequency standards 0, 0, . Of the currently studied 
systems, ^Sq - ^Pi Strontium (Sr) MOTs 6] are partic- 
ularly relevant to fundamental atomic physics since the 
single-photon recoil frequency shift cor is comparable to 
the natural linewidth F and thus ujr directly influences 
both mechanical and thermodynamic trap properties. 
To date, however, many of the rich dynamics for this 
unique system remain experimentally unexplored. 

In this Letter we report a set of novel ^5*0 - ^Pi ^^Sr 
MOT thermal and mechanical dynamics. For laser fre- 
quencies {ujl) tuned below the atomic resonance (cc;^), 
i.e., 27rd=A=ujL — < 0, trap dynamics separate into 
three regimes defined by the relative size of |A|, F, and 
F^;, where F/27r = 7.5 kHz and the power-broadened 
linewidth F^; = F^l + s is determined by the satura- 
tion parameter s = I /Is- Here / {Is = 3 /iW/cm^) 
is the single-beam peak intensity (^5'o - ^Pi saturation 
intensity). Importantly, T ^ uor^ where uor /2t: = 4.7 
kHz. In regime (I), |A| ^ F^; ^ F and semiclassical 
physics dominates. Photon scattering arises predomi- 
nantly from single beams over small, well-defined spatial 
ranges. Gravity also plays an important role as the ratio 
R of the maximum light-induced acceleration vs gravity 
hkV /2mg is only ~ 16, where 27rfi is Planck's constant, 
k is the light wave-vector, m is the ^^Sr mass, and g 
is the gravitational acceleration. Trapped atoms relo- 
cate to vertical positions where magnetic-field-induced 
level shifts compensate \5\ and the resultant radiation 
force balances gravity, leading to ^-independent equilib- 
rium temperatures. In regime (II), |A| < F^;, F^; 
F, a linear restoring force emerges and thermodynam- 
ics reminiscent of ordinary Doppler cooling including 8- 
and s-dependent temperature minima occur, although 
with values globally smaller than standard Doppler the- 
ory predictions. In regime (HI), s approaches unity, the 
photon-recoil-driven impulsive force dominates, and the 
temperature falls below the photon recoil limit {Tr = 



2hjOR/kB = 460 nK, where ks is Boltzmann's constant) 
as predicted by a fully quantum treatment 7]. The fact 
that F ~ ix)R also enables observations of novel J > dy- 
namics, where the ultracold sample divides into momen- 
tum packets whose alignment resembles lattice points on 
3D face-centered cubic crystals. This unique behavior, 
occurring without atomic or excitation coherence, is first 
motivated by an analytic solution to the ID semiclassical 
radiative force equation. Here, we show that for ^ > 0, 
V uor allows direct visualization of "positive feedback" 
acceleration that efficiently bunches the atoms into dis- 
crete, well-defined momentum packets. The experimen- 
tally observed 3D crystal structure is then shown to arise 
naturally from the 3D excitation geometry. In addition, 
we experimentally demonstrate that for fixed 5, 5 de- 
termines the lattice point filling factors, results that are 
confirmed by numerical simulations of the final atomic 
velocity and spatial distributions. More surprisingly, we 
find that R directly influences 5 > {) thermodynamics, 
enabling cooling around a velocity where radiation pres- 
sure and gravity balance. The physics underlying this 
novel behavior is fundamentally the same as regime (I) 
5 < {) cooling, but manifest in a dramatically different 
fashion. 



- ^Pi traps are formed by first pre-cooling ^^Sr in 
a 461 nm - ^Pi MOT with an axial magnetic field 
gradient dB;2/dz (oriented along gravity) of 50 G/cm. 
The atoms are then transferred to 689 nm ^Sq - ^Pi 
MOTs by rapidly lowering dB^/dz to 3 G/cm and ap- 
plying red-detuned broadband frequency-modulated 689 
nm light 6]. Over the next 50 ms, the cloud is com- 
pressed by linearly increasing dB^/dz to 10 G/cm. Sub- 
sequently, highly stabilized, single- frequency 689 nm light 
forms the MOT. The optimal transfer efficiency from ^5*0 
- MOTs to i^'o - ^Pi MOTs is - 30%, giving final 
trap populations of ~ 10''. Typical trap lifetimes and 
spatial densities are ~ 1 s and ~ 5x10^^ cm~^, respec- 
tively. Trap dynamics are monitored either by in-situ or 
time-of-flight (TOF) fluorescence imaging. 



To gain intuitive insight into trap dynamics, we start 



with the semiclassical expression for the force along z, 

hkT 



v,(orVy)(m/s) 
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where s' (> s) signifies contributions from other partic- 
ipating beams and gj = 1.5 (/i) is the ^Pi state Lande 
g-factor (Bohr magneton over h). The force along x (or 
y) is similar to Eq. (1), but without gravity. Figure 
1(a) presents Eq. (1) for dB^/dz = 10 G/cm, s = 
= 248, and a range of S values. The force is displayed 
with respect to position (velocity) in the bottom (upper) 
axis, for Vx or Vy=0 {x or y=0). As S decreases, the 
force makes a clear transition from the Regime (I) | A| 
r^; r isolated form where excitation occurs over two 
separate and well-defined spatial ranges to the Regime 
(II) |A| < r^;, r^; ^ r dispersion-shaped form wherein 
excitation occurs over the entire trap volume and cloud 
dynamics consist of damped harmonic motion 8] . Corre- 
spondingly, as 5 decreases, the initially box-shaped trap 
potential (Fig. 1(b)) becomes progressively more "U"- 
shaped and the trap shifts vertically upward. Finally, 
in regime (III) where F^; approaches F at small 5, sin- 
gle photon recoils dramatically infiuence trap dynamics 
which in turn requires a full quantum treatment J]. 

Changes in the force are dramatically revealed in trap 
mechanical dynamics (see Fig. 1(c)). In the dispersion- 
shaped cooling regime the cloud aspect ratio is ~ 2:1, 
as expected for a typical MOT. Conversely, in the iso- 
lated regime the atoms move freely between "hard wall" 
boundaries. The cloud horizontal width, therefore, is 
largely determined by the separation between horizontal 
force maxima, an effect clearly revealed by the overlaid 
maximum force contours calculated from Eq. (1). More- 
over, since the radiative force is comparable to gravity 
(recall R ~ 16) and the thermal energy is small com- 
pared to the gravitational potential energy, atoms sag to 
the bottom of the trap and the lower cloud boundary 
zq is well defined by the point where the Zeeman shift 
balances S. 

Studying the MOT temperature versus S and s pro- 
vides rich information about trap dynamics. For large 
1^1 and s, corresponding to regime (I), Eq. (1) refiects 
a balance between gravity and the radiative force from 
the upward-propagating beam at 2^0 . Thus trap ther- 
modynamics are determined by a Taylor expansion of 
Eq. (1) around Vz=0 for z=zo. With the atomic posi- 
tion (zq) self-adjusting to follow (5, the damping and mo- 
mentum diffusion coefficients are (5-independent, giving a 
predicted (5-independent temperature of 

T{s) = hTE/{2kB)[0.5R{R - s' / s - l/s)-^^\ (2) 

We have experimentally confirmed this prediction for a 
wide range of 5 10]. The quantity inside the square 
brackets is nearly 2, independent of s for the relevant 
experimental range. Fig. 1(d) displays the temperature 
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FIG. 1: (Color online) (a) Calculated radiative force versus 
position (bottom axis, Vx^Vy^^) and velocity (upper axis, 
x—y—0). (b) Trap potential energy in the z direction, (c) 
In-situ ^5*0 - ^Pi trap images. Dashed lines are calculated 
maximum force contours. For each, s = 248 and dB^^/dz = 
10 G/cm. (d) Temperature vs intensity for 5 — -520 kHz and 
dB^/dz = 10 G/cm. Solid curve: standard Doppler theory; 
Long dashed line: Doppler limit {hV E/'^kB)] Short dashed 
line: single photon recoil limit {2hooR/kB); Filled dots: ex- 
perimental data. 



vs. intensity at a fixed large detuning 6 = -520 kHz, 
showing good agreement (aside from a global scaling fac- 
tor of 2) with the intensity-dependence given by Eq. (2). 
This result arises from the semiclassical nature regime (I) 
cooling for which F^; is the natural energy scale lOj. 

For regime (H), |A| < F^;, F^; > F, Eq. (1) produces a 
linear restoring force resembling ordinary Doppler cool- 
ing. Here, we observe ^3 ^~ 5-dependent tempera- 
ture minima with the minimum and its |^|-location both 
decreasing with s. Such behavior is predicted by Doppler 
theory, with the "Doppler Limit" achieved at |A| = F£;/2. 
However, in order to match the data, the theory curves 
need to be multiplied by a 5-dependent global scaling fac- 
tor (< 1) whose value decreases with s. Moreover, min- 
imum temperatures lie well below the standard Doppler 
limit of hTEf^^B' Notably, this temperature scaling fac- 
tor is not explained by semiclassical Monte-Carlo treat- 
ments of the cooling process. In regime (III), T ^ ur ^ 
kBT/h and the radiative force acquires a single-photon- 
recoil dominated impulsive form. Thus equilibrium ther- 
modynamics can only be adequately described by quan- 
tum theory As shown in Fig. 1(d), the predicted cool- 
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FIG. 2: (a) Underlying momentum space structure for S > 
0. (b) Top view TOF images for in = tv = 25 ms, t^ = 20 
ms. Arrows in the ^ = 60 kHz frame give horizontal trapping 
beam directions, (c) Side view in-situ images for S = 140 kHz 
and Ih — 25 ms. For each, s = 30 and dB^/dz — 0. 



ing limit of half the recoil temperature Tr/2 = huR/kB 
is experimentally reached as s approaches unity. 

Tuning to S > presents another intriguing set of 
cooling and motional dynamics. Here, the cloud divides 
into discrete momentum packets whose alignment mim- 
ics lattice points on a three-dimensional face-centered- 
cubic crystal 0. Figure 2(a) depicts the underlying 
momentum-space structure which, as shown below, oc- 
curs due to highly directional (i.e., minimal heating) 
"positive feedback" acceleration and velocity bunching. 
For the 3D excitation geometry, symmetry dictates that 
cube corners correspond to three-beam processes while 
mid-points between corners and cube face centers arise 
from two- and one-beam processes, respectively. Figure 
2(b) shows a ^-specific sequence of top view (slightly off 
vertical) TOF images for a fixed intensity and atom-light 
interaction time in = 25 ms (ty = 25 ms) in the hori- 
zontal x-y plane (along z-axis) trapping beams, followed 
by a free-flight time tp = 20 ms. All images are taken 
with dBz/dz = although we find qualitatively similar 
behavior for dB^/dz 7^ 0. 

At small values of S (< 60 kHz), the atom cloud ex- 
pands nearly uniformly. As S increases, 3-beam "lattice 
points" appear first, corresponding to the eight cube cor- 
ners in Fig. 2(a). This occurs as the cloud is divided 
into two oppositely moving packets along each of the 



three axes. When S reaches a value around 140 kHz, 
the atom-light interaction becomes sufficiently weak for 
velocities near zero that some atoms remain stationary 
along a given axis. These atoms, however, still inter- 
act with the beams along the two other axes causing the 
2-beam lattice points to appear. This process forms a 
total of 20 divided atom packets with 8, 4, and 8 pack- 
ets present in the top, middle, and bottom layers of the 
cube, respectively. As 5 increases further, some atoms 
are left stationary along two axes, enabling formation of 
the I-beam lattice points, shown as 6 open circles on 
the Fig. 2(a) cube face centers. For S > 180 kHz, the 
atom-light interaction weakens further and the original 
atom cloud reappears. We emphasize that the temper- 
ature associated with each packet in its moving frame 
is actually lower than the ty = tn = atomic cloud. 
This result arises from the velocity bunching and cooling 
mechanisms explained below. 

Only two vertical layers are observed in Fig. 2(b) while 
Fig. 2(a) predicts the creation of three. This apparent 
contradiction is resolved in Fig. 2(c), where the cloud is 
viewed in the x-y plane at 45^ to the x,y axes. In order 
to explore vertical dynamics while maintaining evolution 
in the horizontal plane, in is fixed at 25 ms while ty is 
varied between 6 ms and 25 ms. As before, tp = 20 ms. 
As shown by the images, the lowest two layers in Fig. 2(a) 
are only spatially distinct for short ty, merging together 
for ty = 25 ms. This occurs as gravity accelerates the 
middle layer, which is initially stationary along the z- 
axis, into resonance with the downward propagating laser 
beam. Subsequently the two downward moving layers 
merge. Hence, more (less) intense packets in Fig. 2(b) 
are due to the lowest two (uppermost) cube layers. 

Quantitative insight into S > dynamics can be ob- 
tained from Eq. (1). Recall that for ^ > 0, resonant ab- 
sorption occurs between trapping beams and atoms for 
which k-v > 0. The absorption process thus preferentially 
accelerates rather than decelerates the atoms, leading to 
"positive feedback" in velocity space that terminates at 
a well defined velocity set by s and S. For the unique 
situation where T ^ ujr^ system dynamics rapidly evolve 
toward single-beam interactions. ID dynamics can thus 
be understood by solving Eq. (1) analytically under a 
single-beam approximation. The full 3D evolution then 
follows naturally from the 3D excitation geometry. Fig. 
3(a) illustrates the evolution of the 1-D atomic velocity 
versus interaction time for various 5 at ^ = 100 kHz. Al- 
most independent of the initial velocity {vi), the mean 
value and spread of the final atomic velocity (vj) are set 
by s and 5^ which govern how the acceleration process 
terminates, leading to efficient velocity bunching. Thus 
a S- and 5-dependent number of velocity bunched groups 
are formed. Considering horizontal motion first. Fig 3(b) 
shows, for 5 = 30 and in = 25 ms, Vf versus Vi around 

= for (5 = 100 kHz and S = 160 kHz. In the former 
case, atoms at every Vi are bunched into two groups with 
Vf ^ ± 20 cm/s. In the latter, three groups appear at 
Vf ^ cm/s, 23 cm/s, and -23 cm/s. Similar dynamics 
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FIG. 3: (a) Single beam horizontal acceleration and veloc- 
ity bunching versus time. Numerically calculated two-beam 
final versus initial velocity in the (b) horizontal and (c) ver- 
tical direction. Corresponding spatial distribution in the (d) 
horizontal and (e) vertical direction. 



occur in the vertical direction where gravity now plays an 
important role. Fig. 3(c) shows vj versus Vi for s = 30, 5 
= 140 kHz and ty times relevant to Fig. 2(c). Notably, 
for the upward- moving velocity group, even though S > 
0, atoms experience cooling around a velocity vq where 
gravity balances the radiative force, producing the sharp 
velocity and thus spatial distribution shown in both the 
experiment (Fig. 2(c)) and theory (Fig. 3(c)). Theo- 
retically we find the resultant equilibrium temperature is 



given by Eq. (2), as for the red-detuned case. 

For comparison with Fig. 2, Fig. 3(d) and 3(e) give 
spatial distributions corresponding to the final velocity 
distributions shown in Fig. 3(b) and 3(c), given the 
measured initial cloud temperature. Fig. 3(d) corre- 
sponds to X or y cube axes in Fig. 2(b). Importantly, 
the model correctly reproduces cloud shape asymmetries 
(see, for example, the sharp edge on the top of the up- 
permost layer in Fig. 2(c)), the (5-dependent number of 
packets and the relative packet populations, and the ty- 
dependent number of vertical layers. Predicted final ve- 
locities and packet spacings, however, are ~ 2x larger 
than observed. Measuring the position of the upward 
moving layer in Fig. 2(c) versus ty resolves this discrep- 
ancy. Measured values for Vf are slightly reduced due to 
small stray magnetic field gradients that shift vq as the 
atoms move upward 13], giving an apparent downward 
acceleration. When these effects are taken into account, 
predicted and measured positions agree. Finally, we note 
that ^ > momentum-space crystal formation is a uni- 
versal feature of Doppler limited cooling ^] . For broad 
line cases where T/lor ^ 1, however, creating structures 
similar to Fig. 2(b) requires laser beam diameters on 
the order of tens of centimeters and imaging light with 
hundreds of megahertz bandwidth, making experimental 
observations impractical. 

In summary, we have performed detailed studies of the 
transition from semi-classical to full quantum cooling, re- 
vealing signatures of each regime without ambiguity. Our 
results show, for the first time, that the cooling limit of 
Tr/2 can be reached. More surprisingly, when ^ > 0, the 
cold atom sample divides into well defined momentum 
packets and cooling is achieved around a velocity where 
gravity balances the radiative force. 
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